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Model symulacyjny szacowania niezawodności 
regenerowanych obiektów z wieloparametrową 
oceną stanu technicznego 


Abstract: The paper presents the developed simulation model allowing one to estimate the 
reliability of objects whose technical condition is described by a set of geometrical features. 
Their values change during operation because of the wearing processes of the part under 
consideration. The operational availability (or unavailability) of a technical object is 
determined by the evaluation conducted using a set of criteria referring to its analysed 
features. In addition, the proposed model enables one to take into account the regeneration 
of such objects. 
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Streszczenie: W artykule zaprezentowano opracowany model symulacyjny przeznaczony do 
wyznaczania niezawodności obiektów, których stan techniczny charakteryzowany jest za 
pośrednictwem zespołu cech geometrycznych. Ich wartości ulegają zmianom w czasie 
eksploatacji, na skutek postępującego zużycia. O zdatności (albo niezdatności) obiektu 
technicznego decyduje ocena przeprowadzana z wykorzystaniem zestawu kryteriów 
odnoszących się do jego analizowanych cech. Zbudowany model umożliwia ponadto 
uwzględnienie regeneracji takich obiektów. 
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1. Introduction 


Simulation methods, despite their disadvantages, are used in research of 
various systems: physical, military, social, energetic, political and economic. The 
complexity of these systems, and in particular the difficulties in obtaining solutions 
using analytical methods, necessitate the use of simulation methods. Such a wide 
range of applications results mainly from the possibility of using a number of types 
of simulations and their adaptation to various models of the reality being examined 
[2, 7, 8, 14]. 

Due to its advantages, digital simulation is also used in reliability analyses. 
The use of simulations for estimating the reliability of renewable systems with 
structural redundancy in the form of additional elements is presented in [5]. The 
influence of the number of redundant elements on the reliability of the entire system 
was demonstrated in cases when the system consists of a large number of basic 
elements. The study [9] presents the use of simulation to estimate the reliability of 
a complex technical system with a structure of type k out of n, in which the elements 
are renewed after failure and can be renewed preventively if they have been 
working without failure during a fixed period of time. The introduced changes of 
time till preventive renewal of the elements cause changes in the reliability of the 
system in the considered period. The results obtained in simulation experiments 
allow to determine such time till preventive renewal of elements, so that in the 
analysed time the system can be characterised by the assumed level of reliability. 
Stochastic simulation is also used in estimating the reliability of systems having 
network structures in which there is a problem of determining the reliability 
structure (cf. [6]). The study [10] presents the possibility of estimation of reliability 
of modelled networks using simulation, graph theory and the factorization 
algorithm. Obtained results can be used in the selection of network structures at the 
stage of their design, modernisation and extending. 

The paper presents a developed simulation calculation model for estimating 
the reliability of regenerated objects with a multi-parameter assessment of their 
technical condition. The correct functioning of this model has been verified on the 
example of an object, which is the wheel of a railway vehicle. 


2. Multistate regenerated technical objects 


It has been assumed that the reliability state of a technical object, which may 
be, for example, an analysed in the article rail vehicle wheel, depends on its 
technical condition, and the object itself may be available (when it may work 
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properly) or unavailable (otherwise). At the same time, the technical condition of 
the object described with n significant features due to its proper operation is 
characterised by a feature vector (cf. [12]): 


X(t, w) = [xw (t), X2w(t), Ese) Xn,w OJ 


where: 
Xiw(t) — value of feature i of the object operating in the time period (0, t) in 
conditions w, 
t — moment in the history of operation and maintenance process in which the 
technical condition of the object is described, 
— number of considered features of the object, 
w — operation and maintenance conditions in the time period (0, t). 


a 


Due to the conditions of the operation and maintenance practice, it is important 
to determine the relationship between the technical condition of the object and its 
reliability state. It was assumed that the whole object remains available if and only 
if the value of each of its distinguished features does not exceed the accepted 
(required) limit values [12, 13]: 

— for one feature with the value of c(ż) at time t: 


OT a(t) = [(c(t) = Cmin) A (c(t) < Cmax)I, 


— for n features with the values c;(/) at times ¢, i = 1, 2, ..., n: 


OTza(t) © NiCCmin < Ci (t) < Cmax). 


Therefore, it may be said that in the operation and maintenance process, the 
examined technical object may become unavailable at time ¢ (OTnza(f)) due to 
exceeding by the value of at least one of its considered features of the appropriate 
limit value (minimum or maximum). This is correct on the assumption that as a 
result of the course of the object's operation and maintenance process, the values 
of its features may change in a continuous or stepwise manner. This is a typical 
undesirable consequence of the progressive degradation processes (e.g. wear, 
fatigue, corrosion, erosion, aging) accompanying its use. 

With such a description of the technical condition of a technical object and the 
formulation of the condition of its availability, the desired changes in the technical 
condition of the object caused by its regeneration can be interpreted as changes in 
the value of its certain features. They occur by increasing the absolute value of the 
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difference between the value of i feature and the relevant limit value and lead to 
extending the durability (cf. [12, 13]): 


[lcir (t) — Ciminl > [Ciz (t) — Ciminl], 


or 


[lcir Ct) m" Cimaxl > lciz (t) = Cimaxl], 


where: 
r,z — relate respectively to the values of the features after and before the 
regeneration of the object operating in the time period (0, t). 


The proper description (one out of the two mentioned cases) may be applied 
depending on whether the value of its i feature decreases or increases as a result of 
using the object. 

This problem was solved in the article using the developed simulation 
program. 


3. Simulation model for reliability estimation 


The reliability estimations obtained with the use of computer simulation allow 
to support decision-making processes in the range of rational choices regarding: 
system’s structure, type of redundancies used, maintenance system, number of 
redundant objects or strategies of inspections and preventive renewal of objects 
used. They are already made at the design stage of the technical system. These 
issues were analysed and presented in many publications [1, 4, 5, 15]. 

This article presents a simulation method for evaluation of the reliability of 
objects characterised by the fact that their wear process is proceeding in several 
different ways and each of them causes changes in the value of specific features. 
Each of the features has an accepted limit value, the achievement of which is 
tantamount to the failure of the object due to this feature. In accordance with the 
activities applied in the operational practice of these objects, it has been assumed 
that it is possible to regenerate the object several times during its operation. 
Conducted regeneration process causes a step change in the value of specific 
features — to the initial state (like new), at the expense of increasing the degree of 
wear of the object characterised by another feature, which determines the 
possibility of further regeneration of the whole object. The values of this feature 
are also subject to changes due to wear during the operation, regardless of 
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regeneration. Therefore, the technical condition of the object after regeneration is 
not equivalent to the technical condition of the object after complete renewal (like 
new), although it is assumed that the values of some of its features correspond to 
the situation of the object fully renewed. The overall reliability structure of such an 
object is shown in fig. 1. 


Regenerated elements Ep Non-regenerated elements Ey 


a 


Fig. 1. Overall reliability structure of analysed object (reflecting the dependence of 
object’s reliability state on its selected features) 


On the basis of information about the real operation process of such objects, a 
simulation model was developed that reflects the random process of the object's 
wear in view of each of the distinguished features and includes carrying out the 
regeneration of the object until its total durability expires. 

The iteration algorithm is shown in fig. 2. 

The input to the algorithm is: reliability structure of a system, specified number 
of regenerations, reliability functions of system elements, time horizon and time 
step. 

To implement the task, stochastic simulation with a constant time step has been 
used, in which random moments (in accordance with the probability distribution) 
of the object’s failures are determined due to exceeding the value of a given feature. 
If the failure is the result of exceeding the value of the factor that determines 
conducting regeneration process and their number has not exceeded the limit, 
regeneration restores the object's availability and initial state of the relevant 
features, and the iteration of the simulation continues. If the failure is the result of 
exceeding the value of the feature not decisive for regeneration, then the object 
becomes unavailable and the iteration is interrupted. A single iteration of the 
simulation ends with the object’s failure or is successful if the object does not 
become unavailable before the time horizon expires. 

The simulation experiment includes the implementation of many iterations. In 
practice, their number is a minimum of 10,000. 


11 


Stanistaw Mtynarski, Robert Pilch, Maksymilian Smolnik, Jan Szybka, Grzegorz Wigzania 


Conducting experiments according to the developed simulation allows 
estimating the reliability of the considered objects in the assumed time horizon of 
operation and maintenance, which can be used in planning the times of their 
inspections, exchanges and the necessary number of spare parts. 


Input data: 
1. Structure of the system, 

Ep - regenerated elements, 

En - non regenerated elements 
Max. number of regenerations: Np 
CDF of elements 

Time horizon: Tp 
Time step: dt 


PUUN 


tu> Th 
or 
ig > Nę 1. Computing values: Pug gy and Pug p, 
or 2. Drawing RND: pj and RNDE yi 
for each Ex, 3. Increasing age of elements te gy, te ni 
RNDE w < PUE nx and the system ty by adding dt 


YES 


Save age of 
the system ty 


NO for each Eg; 


RNDE pj < PUE pj 


Fig. 2. Algorithm of a single iteration of the developed simulation procedure 


4. Physical aspects of rail vehicles wheels wear 
The estimation of reliability indexes of objects whose parts are subject to 


degradation at different paces with each other is of great importance in the analysis 
of reliability and durability of real technical objects. A good example of this type 
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of objects are wheels of railway vehicles, on the example of which the analysed 
problem is presented below. 

The safety of rail vehicle operation is primarily dependent on the availability 
of wheelsets. The intensity of wear and the failure rate of sudden failures resulting 
from random factors damaging vehicle wheels is of paramount importance in the 
use of rolling stock [3, 9]. 

Ensuring the safety of the operation of rail vehicles depends on the 
preservation of a wheel profile set out by the standard. The wheel profile is a 
specially shaped curvature of the outer surface of the wheel rim, which consists of 
two parts: a tread surface, wheel flange. Meeting the requirements of the standard 
is accomplished by maintaining specific boundary dimensions of both the tread 
surface and the wheel flange. A view of an exemplary wheel profile [16] is shown 
in fig. 3 [17]. 

Wheels of rail vehicle wheelsets still constitute a significant maintenance 
problem and in use they may reach limit values of wear and are subjected to damage 
and require frequent profile renewal. Analysis of the characteristics provides data 
on the behavior of these elements in the conditions of their operation. 

During the operation of vehicles there is wear of wheelsets on the wheel tread, 
and in particular the area in the region of the flange constituting the transition zone 
to the traed. The observation of the operation and maintenance process shows that 
the vehicles in the service cycle, in accordance with the maintenance system 
documentation, during the mileage of 500,000 km undergo two profile regeneration 
treatments (so-called full reprofiling) with parameters defined for the S1002 profile 
according to the PN-EN 13715 standard [16]. 

The comparison of factors affecting the durability of wheelsets in the aspect of 
wheel flange wear and wheel wear on the running thread in the case under 
consideration indicates that wear of the flange (‘undercut’ or sharp flange) is a 
decisive factor in the need to reprofile wheels. The running thread wear is not a 
factor determining the durability of a wheelset, despite the fact that the limit size of 
the running thread determines the withdrawal of the wheelset from operation. The 
dominant influence on reducing the value of the ‘rim diameter’ feature is the fact 
that the edge is reprofiled. 

The analysed results of operational tests of wear showed that ‘under-cutting’ 
of the wheel flange determines the need for regeneration of the majority of 
wheelsets. The renewal of the tread profile and dimensions of the flange depends 
on the dimensions of the tread and, in order to maintain the required dimensions of 
the wheels — of the wheelset and the dimensions of all wheels of the vehicle trolley 
determines the durability of the wheelsets. Taking into account the results of the 
tests, it is also important to accurately plan the renewal operations of the rail vehicle 
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trolleys and forecast the durability of the wheel rims. The obtained results of the 
simulation of durability and reliability of rail vehicle wheels with the proposed 
method can be used to optimise the maintenance process and improve the 
economics of railway vehicle operation and other technical objects with 
multiparameter evaluation of their technical condition. 


Fig. 3. A view of an exemplary wheel profile [17] — marking of measurable features of the 
rolling profile of the wheel: G — tyre thickness, Go — flange thickness, Wo — flange 
height, qr — flange slope, b — width of the tyre or the monoblock wheel rim, Az — 
distance between tyre back surfaces or rims of monoblock wheels in the wheelsets 
without load, D — rim diameter, D,— wheel centre diameter, Ez — distance between 
flange profiles 


5. Exemplary calculation 


In the analysed case, based on the data collected during the operation of 
selected rail vehicles, the distributions of the probability of mileage to reaching the 
limit state were determined. The object under consideration is the wheel of a rail 
vehicle. Selected features describing the wear (and also the technical condition) of 
the object are: wheel diameter (D), flange thickness (GO) and wheel profile outline 
(P). The obtained probability density functions (f(s)) of mileage to reaching the 
limit state for each of the analysed features are shown in fig. 4. 

The diameter and outline of the profile are described using Weibull distribution 
with a density function in the form of: 
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and the outline of the profile using normal distribution. The values of the 
distribution parameters are given in tab. 1. 


Table 1 
The values of the distribution parameters 
Wheel feature Profile outline P Flange thickness Go Wheel diameter D 
Facies normal: Weibull: Weibull: 
SŁAWKA m= 215000 km v=5 v=3,4 
its parameters 


6 = 60000 km B= 158000 km B = 429000 km 


Fig. 4. Graphs of probability density functions (f(s)) of mileage to reaching the limit state: 
fD(s) — wheel diameter, fP(s) — profile outline, fGo(s) — flange thickness 


The wheel reliability structure, determined by a set of listed features, is shown 
in fig. 5. 
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Regenerated elements Ep 


Non-regenerated 
elements 
En 
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Fig. 5. Considered reliability structure of a rail vehicle wheel 


On the basis of the adopted model of the system (fig. 5) and the algorithm of 
the simulation procedure (fig. 2) a computer program was developed to simulate 
system’s operation. The simulation parameters were: 


— time step dt equal to 5000 km, 
— time horizon TH equal to 600000 km, 
— number of iterations n equal to 10000. 


The reliability estimations obtained from the conducted simulations are shown 


in fig. 6. 
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Fig. 6. Obtained courses of reliability function of a rail vehicle wheel: R(S) — reliability of 
the wheel when regeneration is not carried out, R1(s) with one and R2(s) with two 


regenerations 
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The analysis of the obtained mileage indicates a significant effect of 
regeneration on the reliability of operation and durability of the considered vehicle 
wheel. Culmination of failures by mileage of about 150000 km (fig. 4) now does 
not directly cause a significant change in the reliability of the entire object after it 
reaches this mileage (s). When only one regeneration is performed, the mileage is 
almost doubled to reach a similar level of reliability (R2(s)), and with double 
regeneration almost three times (R3(s)). However, a non-regenerated feature, which 
is the diameter of the wheel, has a significant impact on reliability. 


6. Conclusions 


In the field of analysis and estimation of reliability and readiness of technical 
systems, the most common tool applied is a stochastic computer simulation, an 
example of which is presented in the article. This simulation makes it possible to 
obtain reliability estimations in much shorter time than would be necessary to 
evaluate reliability on the basis of observations of real systems and with 
significantly lower costs. Due to the possibility of repeating the experiment in 
unchanged conditions, it allows to obtain large numbers of parameter values (which 
can be referred to as values of considered features) required to evaluation reliability, 
and difficult to observe in practice, especially when the analysed systems are 
unique or present in a small number. In addition, the values of features that are 
important in the evaluation of reliability, and which are random variables of 
specific probability distributions (such as working times till failure, duration of 
renewal, etc.) in a computer simulation experiment may easily be obtained by using 
random number generators. 

The simulation method presented in the paper may be successfully used to 
predict the durability of rail vehicle wheels. In the developed simulation model, the 
influence of regeneration activities (wheel reprofilation) on extending the working 
time of wheelsets of rail vehicles is evident. It is also possible to use it for planning 
and synchronising regenerative activities of individual wheels functioning in wheel 
sets. 


This work was financed by AGH University of Science and Technology, Faculty of 
Mechanical Engineering and Robotics, research program No. 11.11.130.174. 
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MODEL SYMULACYJNY SZACOWANIA 
NIEZAWODNOSCI REGENEROWANYCH OBIEKTOW 
Z WIELOPARAMETROWĄ OCENĄ STANU 
TECHNICZNEGO 


1. Wprowadzenie 


Metody symulacyjne, mimo wad, znajdują zastosowanie w badaniach różnych 
systemów: fizycznych, wojskowych, społecznych, energetycznych, politycznych, 
a także ekonomicznych. Złożoność tych systemów, a w szczególności trudności 
w uzyskaniu poszukiwanych rozwiązań w drodze analitycznej, powoduje 
konieczność stosowania metod symulacyjnych. Tak szeroki obszar zastosowań 
wynika przede wszystkim z możliwości wykorzystywania wielu typów symulacji 
i ich adaptacji do różnych modeli badanej rzeczywistości [2, 7, 8, 14]. 

Ze względu na swoje zalety symulacja cyfrowa jest również wykorzystywana 
w analizach niezawodności. Zastosowanie symulacji do szacowania niezawodności 
odnawialnych układów posiadających nadmiar strukturalny w postaci 
dodatkowych elementów przedstawiono w [5]. Wykazano wpływ liczby 
elementów rezerwowych na niezawodność całego układu, w przypadkach gdy 
układ składa się z dużej liczby elementów podstawowych. W opracowaniu [9] 
zaprezentowano wykorzystanie symulacji do wyznaczania niezawodności 
złożonego układu technicznego o strukturze typu k z n, w którym elementy są 
odnawiane po uszkodzeniu oraz mogą być odnawiane prewencyjnie, jeśli 
bezawaryjnie przepracowały ustalony czas. Wprowadzane zmiany okresów 
odnowy prewencyjnej elementów powodują zmiany niezawodności układu 
w przyjętym okresie. Uzyskane w eksperymentach symulacyjnych wyniki 
umożliwiają wyznaczenie takich okresów odnów prewencyjnych elementów, aby 
w analizowanym czasie układ utrzymał założony poziom niezawodności. 
Symulacja stochastyczna jest również wykorzystywana w  szacowaniu 
niezawodności układów posiadających struktury sieciowe, w których występuje 
problem określenia struktury niezawodnościowej (por. [6]). W opracowaniu [10] 
przedstawiono możliwość symulacyjnego szacowania niezawodności sieci 
modelowanych z wykorzystaniem teorii grafów i algorytmu faktoryzacji. 
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Uzyskiwane wyniki mogą być wykorzystane w doborze struktur sieci na etapie ich 
projektowania, modernizacji i rozbudowy. 

W artykule pokazano opracowany symulacyjny model obliczeniowy 
przeznaczony do szacowania niezawodności regenerowanych obiektów 
z wieloparametrową oceną ich stanu technicznego. Poprawność funkcjonowania 
tego modelu zweryfikowano na przykładzie obiektu, jakim jest koło pojazdu 
szynowego. 


2. Wielostanowe regenerowane obiekty techniczne 


Przyjęto założenie, że stan niezawodnościowy obiektu technicznego, którym 
może być np. analizowane w artykule koło pojazdu szynowego, jest zależny od 
jego stanu technicznego, a sam obiekt może być zdatny (gdy może pracować 
poprawnie) albo niezdatny (w przeciwnym przypadku). Jednocześnie, stan 
techniczny obiektu opisywanego n cechami istotnymi z uwagi na jego prawidłową 
pracę, charakteryzuje wektor cech (por. [12]): 


X(t, w) = iw (t), X2w(t), ese) Xn,w OJ 


gdzie: 

Xiw(t) — wartość i-tej cechy obiektu eksploatowanego w czasie (0, t) w warun- 
kach w, 

t — chwila w historii eksploatacji, w której opisywany jest stan techniczny 
obiektu, 

n — liczba rozpatrywanych cech obiektu, 

w — warunki eksploatacji obiektu w czasie (0, t). 


Ze względu na uwarunkowania praktyki eksploatacyjnej, istotne jest zatem 
określenie zależności pomiędzy stanem technicznym obiektu i jego stanem 
niezawodnościowym. Założono, że cały obiekt pozostaje zdatny wtedy i tylko 
wtedy, gdy wartość każdej z jego wyróżnianych cech nie przekracza przyjętych 
(wymaganych) granicznych wartości [12, 13]: 

— dla jednej cechy o wartości c(f) w chwili t: 


OTza(t) S [(C(t) 2 min) A (C(O) £ Cmax)], 


— dla n cech o wartościach c:( w chwili t, i= 1, 2, ..., n: 
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OTza(t) © NiCCmin < Ci (t) < Cmax). 


Tym samym można powiedzieć, że w procesie eksploatacji badany obiekt 
techniczny może stać sie niezdatny w chwili £ (OTnza(t)) z powodu przekroczenia 
przez wartość choćby jednej, dowolnej z rozpatrywanych jego cech odpowiedniej 
dla niej wartości granicznej (minimalnej albo maksymalnej). Jest to słuszne przy 
założeniu, że na skutek przebiegu procesu eksploatacji obiektu, wartości jego cech 
mogą zmieniać się w sposób ciągły lub skokowy. Jest to typowym niepożądanym 
następstwem postępujących procesów jego degradacji (np. zużycia, zmęczenia, 
korozji, erozji, starzenia), towarzyszących jego wykorzystywaniu. 

Przy tak dokonanym opisie stanu technicznego obiektu technicznego oraz 
sformułowaniu warunku jego zdatności, pożądane zmiany stanu technicznego 
obiektu spowodowane jego regeneracją mogą być interpretowane jako zmiany 
wartości jego pewnych cech. Zachodzą one, zwiększając wartość bezwzględną 
różnicy wartości i-tej cechy i odpowiedniej wartości granicznej oraz prowadzą do 
wydłużenia czasu eksploatacji (por. [12, 13]): 


[lcir Ct) — Cimin| > lciz (t) — Cimin l], 


albo 


[lcir Ct) r Cimaxl > lciz (t) = Cimaxl], 


gdzie: 
r,z — dotyczą odpowiednio wartości cech po i przed regeneracją obiektu 
eksploatowanego w czasie (0, t). 

Właściwy opis spośród dwóch wymienionych przypadków znajduje 
zastosowanie zależnie od tego, czy na skutek wykorzystywania obiektu wartość 
jego i-tej cechy maleje, czy wzrasta. 

Problem ten postarano się rozwiązać w artykule, stosując autorski program 
symulacyjny. 


3. Symulacyjny model szacowania niezawodności 
Oszacowania niezawodności uzyskiwane z zastosowaniem symulacji 
komputerowej pozwalają wspomagać procesy decyzyjne w zakresie racjonalnych 


wyborów dotyczących: struktury układu, rodzaju stosowanych nadmiarów, 
systemu obsługiwania, liczby obiektów rezerwowych, czy też strategii przeglądów 
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i odnów prewencyjnych obiektów. Dokonywane są one już na etapie projektowania 
układu technicznego. Zagadnienia te analizowano i prezentowano w wielu 
publikacjach [1, 4, 5, 15]. 

W niniejszym artykule przedstawiono opracowaną symulacyjną metodę oceny 
niezawodności obiektów charakteryzujących się tym, że ich zużycie przebiega na 
kilka różnych sposobów i każdy z nich powoduje zmiany wartości określonych 
cech. Każda z cech posiada przyjętą wartość graniczną, której osiągnięcie jest 
równoznaczne z uszkodzeniem obiektu ze względu na tę cechę. Zgodnie 
z działaniami stosowanymi w praktyce eksploatacyjnej tych obiektów przyjęto, że 
istnieje możliwość kilkukrotnej regeneracji obiektu w okresie jego eksploatacji. 
Przeprowadzana regeneracja powoduje skokową zmianę wartości określonych 
cech — do stanu początkowego (jak nowe), kosztem zwiększenia stopnia zużycia 
obiektu charakteryzowanego inną cechą, która warunkuje możliwość kolejnych 
regeneracji całego obiektu. Wartości tej cechy ulegają również zmianom na skutek 
zużycia postępującego w trakcie eksploatacji niezależnie od wykonywanych 
regeneracji. Zatem stan techniczny obiektu po wykonaniu regeneracji nie jest 
równoważny stanowi technicznemu obiektu po całkowitej odnowie (jak nowy), 
mimo to przyjmuje się, że wartości pewnych jego cech odpowiadają sytuacji 
obiektu w pełni odnowionego. Ogólną strukturę niezawodnościową takiego obiektu 
pokazano na rys. 1. 


Elementy regenerowane Ep Elementy nieregenerowane Ey 
r Eri ta Eg a Er 7 Ena Kex Enk S Enm 
n a 
n n 
" n 
" a 
Ww 1" 
" ti 
yk. 11 
i Ye 11 
my. i h 
tt 
' zj Eri Ej Er T I 
1 1 
1 E 1 
1 1 
1 1 
U 1 
U 
L ER bas ERj Pasz Eni a! 


Rys. 1. Ogólna struktura niezawodnościowa analizowanego obiektu (uzalezniajaca stan 
niezawodnościowy obiektu od jego wyróżnianych cech) 


Na podstawie informacji o rzeczywistym procesie eksploatacji takich 
obiektów opracowano model symulacyjny, który odzwierciedla losowy proces 
zużywania się obiektu ze względu na każdą z wyróżnionych cech i uwzględnia 
przeprowadzanie regeneracji obiektu aż do osiągnięcia całkowitego jego zużycia. 

Algorytm iteracji przedstawiono na rys. 2. 
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Dane wejsciowe: 
1. Struktura układu, 

Erę - elementy regenerowane, 

Ey - elementy nieregenerowane 
2. Graniczna liczba regeneracji Nz 
3. Funkcje niezawodności elementów 
4 
5 


. Horyzont czasowy: Ty 
. Krok czasowy dt 


tu > TH 
lub 
ip > Nę 
lub 
dla dowolnego Ex, 
RNDe nk < PUE nk 


1. Wyznaczenie wartości Pug pj i PUE ni 
2. Wylosowanie liczb RNDz pj Í RNDE ji 
3. Zwiększenie wieku elementów 

te rj te ni i układu ty o wartość dt 


TAK 


Y 
Zapisanie wieku 
układu ty 


Rys. 2. Algorytm pojedynczej iteracji opracowanej procedury symulacyjnej 


dla dowolnego Ep, 
RND: rj < PUE pj 


Danymi wejściowymi do algorytmu są: struktura niezawodnościowa układu, 
graniczna liczba regeneracji, funkcje niezawodności elementów układu, horyzont 
czasowy oraz krok czasowy. 

Do realizacji zadania wykorzystano symulację stochastyczną, o stałym kroku 
czasowym, w której wyznaczane są losowe (zgodnie z rozkładem prawdopodobień- 
stwa) chwile uszkodzeń obiektu ze względu na przekroczenie wartości określonej 
cechy. Jeśli uszkodzenie jest efektem przekroczenia wartości cechy decydującej 
o wykonaniu regeneracji i ich liczba nie została wyczerpana, następuje regeneracja 
przywracająca stan zdatności obiektu oraz stan początkowy odpowiednich cech, 
aiteracja symulacji trwa dalej. Jeśli uszkodzenie jest efektem przekroczenia 
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wartości cechy niedecydującej o regeneracji, to następuje uszkodzenie obiektu 
iprzerwanie iteracji. Pojedyncza iteracja symulacji kończy się uszkodzeniem 
obiektu lub kończy się sukcesem, jeśli przed upływem horyzontu czasowego 
uszkodzenie obiektu nie nastąpi. 

Eksperyment symulacyjny obejmuje wykonanie wielu iteracji. W praktyce ich 
liczba wynosi minimum 10 000. 

Przeprowadzenie eksperymentów według opracowanej symulacji umożliwia 
wyznaczenie niezawodności rozważanych obiektów w przyjętym horyzoncie 
czasowym eksploatacji, co może być wykorzystywane w planowaniu czasów 
dokonywania ich przeglądów, wymian oraz niezbędnej liczby obiektów 
wymiennych. 


4. Fizykalne aspekty zużycia kół pojazdów szynowych 


Duże znaczenie w analizie niezawodności i trwałości rzeczywistych obiektów 
technicznych ma ocena wskaźników niezawodności obiektów, których elementy 
podlegają degradacji w różnym tempie względem siebie. Dobrym przykładem 
obiektów tego rodzaju są koła pojazdów szynowych, na których przykładzie 
poniżej przedstawiono analizowany problem. 

Bezpieczeństwo eksploatacji pojazdów szynowych przede wszystkim 
uzależnione jest od zdatności zestawów kołowych. Intensywność zużycia oraz 
intensywność uszkodzeń nagłych, wynikających z losowych czynników uszkadza- 
jących koła pojazdów, ma nadrzędne znaczenie w użytkowaniu taboru [3, 9]. 

Zapewnienie bezpieczeństwa eksploatacji pojazdów szynowych uzależnione 
jest od zachowania określonego normą profilu zarysu kół zestawów kołowych. 
Profil toczny koła stanowi specjalnie ukształtowana krzywizna zewnętrznej 
powierzchni wieńca koła, która składa się z dwóch części: powierzchni tocznej 
iobrzeża koła. Wymogi normy spełniane są przez zachowanie określonych 
wymiarów granicznych obydwu powierzchni, tj. tocznej i obrzeża koła. Widok 
przykładowego profilu koła [16] pokazano na rys. 3 [17]. 

Koła zestawów kołowych pojazdów szynowych ciągle stanowią istotny 
problem eksploatacyjny i w trakcie użytkowania mogą osiągnąć graniczne wartości 
zużycia oraz podlegać uszkodzeniom i wymagają zabiegów częstej odnowy zarysu. 
Analiza charakterystyk dostarcza informacji o zachowaniu się tych elementów 
w warunkach panujących podczas ich pracy. 

W trakcie eksploatacji pojazdów następuje zużycie zestawów kołowych na 
powierzchni tocznej obręczy, a w szczególności w rejonie obrzeża, w strefie 
przejścia w powierzchnię toczną. Z, obserwacji procesu eksploatacji wynika, że 
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pojazdy w cyklu miedzynaprawczym zgodnie z DSU (dokumentacja systemu 
utrzymania) w trakcie pokonywania przebiegu 500 000 km maja wykonywane dwa 
zabiegi odtworzenia profilu (tzw. pełnej reprofilacji) o parametrach określonych 
dla profilu S1002 wg normy PN-EN 13715 [16]. 


Rys. 3. Widok fragmentu koła [17] — oznaczenie cech mierzalnych profilu tocznego koła: 
G — grubość obręczy, Go — grubość obrzeża, Wo — wysokość obrzeża, 
qR — stromość obrzeża, b — szerokość obręczy lub wieńca koła bezobręczowego, 
Az — odległość między wewnętrznymi powierzchniami obręczy lub wieńcami kół 
bezobręczowych w zestawach kołowych bez obciążenia, D — średnica okręgu 
tocznego, Dı- średnica koła bosego, E, — odległość między zarysami obrzeży 


Z porównania czynników wpływających na trwałość zestawów kołowych 
pojazdu w aspekcie zużycia obrzeża koła oraz zużycia koła na kręgu tocznym 
w rozpatrywanym przypadku wynika, że zużycie obrzeża „podcięcie” jest 
czynnikiem decydującym o konieczności reprofilacji kół. Zużycie na kręgu 
tocznym nie jest w tym przypadku wielkością decydującą o trwałości zestawu 
kołowego, pomimo że graniczne wymiary średnicy okręgu tocznego decydują 
o wycofaniu zestawu kołowego z eksploatacji. Dominujący wpływ na zmniejszenie 
wartości cechy „średnica na okręgu tocznym” ma fakt wykonywania reprofilacji 
obrzeża. 

Przeanalizowane wyniki badań zużycia eksploatacyjnego wykazały, że 
podcięcie obrzeża obręczy determinuje zużycie większości zestawów kołowych. 
Odnowa profilu i wymiarów obrzeża jest zależna od wymiarów powierzchni 
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tocznej i ze względu na zachowanie wymaganych wymiarów kół — w zestawie 
kołowym oraz wymiarów wszystkich kół wózka pojazdu determinuje trwałość 
zestawów kołowych. Uwzględniając wyniki badań, istotna jest również możliwość 
dokładnego planowania zabiegów naprawczych wózków pojazdu szynowego oraz 
prognozowania trwałości obręczy kół. Uzyskane wyniki symulacji trwałości 
iniezawodności kół pojazdów szynowych proponowaną metodą mogą być 
wykorzystane do optymalizacji procesu obsługi oraz poprawy ekonomiki 
eksploatacji pojazdów szynowych oraz innych obiektów z wieloparametrową 
oceną stanu technicznego. 


5. Przykład obliczeniowy 


W analizowanym przypadku, bazując na danych zgromadzonych w czasie 
eksploatacji wybranych pojazdów szynowych, określono rozkłady prawdopodo- 
bieństwa przebiegu do osiągnięcia stanu granicznego. Rozpatrywanym obiektem 
jest koło pojazdu szynowego. Wyróżnionymi cechami opisującymi zużycie 
(a zarazem stan techniczny) obiektu są: średnica koła (D), grubość obrzeża (Go) 
oraz zarys profilu (P). Uzyskane funkcje gęstości prawdopodobieństwa przebiegu 
(f(s)) do osiągnięcia stanu granicznego dla każdej z analizowanych cech 
przedstawiono na rys. 4. 

Średnicę oraz zarys profilu opisuje rozkład Weibulla z funkcją gęstości 
w postaci: 


ANE sc SU 
194) eG) 
p B 
a zarys profilu rozkład normalny. Wartości parametrów rozkładów zamieszczono 
w tab. 1. 


Tabela 1 


Parametry rozkładów prawdopodobieństwa 


Cecha koła Zarys profilu P Grubość obrzeża Go Średnica koła D 


normalny: Weibulla: Weibulla: 
PARA m= 215 000 km v=5 v=3,4 
eRe ane o = 60 000 km B = 158 000 km B = 429 000 km 
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0,014 


0,012 


0,01 


Rys. 4. Wykresy funkcji gęstości prawdopodobieństwa przebiegu do osiągnięcia stanu 
granicznego: fD(s) — średnica koła, fP(s) — geometria profilu, fGo(s) — grubość 
obrzeża 


Strukturę niezawodnościową koła, wyznaczaną przez zespół wymienionych 
cech, pokazano na rys. 5. 
Element 


nieregenerowany 
Elementy regenerowane Eg Ey 


>———Ą 


Rys. 5. Przyjęta struktura niezawodnościowa koła pojazdu szynowego 


Na podstawie przyjętego modelu układu (rys. 5) oraz algorytmu procedury 
symulacyjnej (rys. 2) opracowano program komputerowy pozwalający na 
symulację jego pracy. Parametry symulacji wynosiły: 

— krok czasowy dt równy 5 tys. km, 

— horyzont czasowy Ty równy 600 tys. km, 


2 


N 


Stanistaw Mtynarski, Robert Pilch, Maksymilian Smolnik, Jan Szybka, Grzegorz Wigzania 


— liczba iteracji n równa 10 000. 


Oszacowania niezawodności uzyskane na podstawie przeprowadzonych 
symulacji przedstawiono na rys. 6. 


NiezawodnośćR(s) 


0 100 200 300 400 500 600 700 
s [103 km] 


Rys. 6. Uzyskane przebiegi funkcji niezawodności koła pojazdu szynowego: R(s) — 
niezawodność koła, gdy nie jest wykonywana regeneracja, R1(s) z jedną, a R2(s) 
z dwiema regeneracjami 


Analiza uzyskanego przebiegu wskazuje na znaczący wpływ regeneracji na 
niezawodność funkcjonowania i trwałość rozpatrywanego koła pojazdu szyno- 
wego. Kulminacja uszkodzeń przypadająca na przebieg ok. 150 tys. km (rys. 4) nie 
przekłada się teraz bezpośrednio na istotną zmianę niezawodności całego obiektu 
po osiągnięciu przez niego tej wartości przebiegu (s). Przy wykonywaniu tylko 
jednej regeneracji obserwuje się prawie dwukrotny wzrost przebiegu do osiągnięcia 
zbliżonego poziomu niezawodności (R2(S)), a przy dwukrotnej regeneracji prawie 
trzykrotny (R3(s)). Istotny wpływ na niezawodność wykazuje natomiast cecha 
nieregenerowana, którą jest średnica koła. 


6. Podsumowanie 
W zakresie analiz i szacowania niezawodności oraz gotowości układów 


technicznych najczęściej wykorzystywana jest stochastyczna symulacja 
komputerowa, której przykład zastosowania przedstawiono w artykule. Symulacja 
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ta umożliwia uzyskanie oszacowań niezawodności w czasie znacznie krótszym niż 
byłby konieczny do oceny niezawodności na podstawie obserwacji układów 
rzeczywistych i przy zdecydowanie niższych kosztach. Dzięki możliwości 
wielokrotnego powtarzania eksperymentu w niezmienionych warunkach, pozwala 
uzyskiwać duże liczności wartości parametrów (które mogą być określane jako 
wartości rozpatrywanych cech) wymaganych do oceny niezawodności, a trudnych 
do zaobserwowania w praktyce, szczególnie gdy analizowane układy eksploato- 
wane są jako unikatowe lub występują w niewielkiej liczbie. Ponadto wartości 
cech, istotnych w ocenie niezawodności, a będących zmiennymi losowymi 
o określonych rozkładach prawdopodobieństwa (jak np. czasy pracy do uszko- 
dzenia, trwania odnowy itp.), w komputerowym eksperymencie symulacyjnym 
można łatwo uzyskiwać poprzez wykorzystanie generatorów liczb losowych. 

Przedstawiona w opracowaniu metoda symulacji może być z powodzeniem 
zastosowana do prognozowania trwałości kół pojazdów szynowych. W opracowa- 
nym modelu symulacyjnym widoczny jest wpływ zabiegów regeneracyjnych 
(odtwarzania profilu koła) na wydłużenie czasu pracy zestawów kołowych 
pojazdów szynowych. Możliwe jest również wykorzystanie jej do planowania 
i synchronizowania zabiegów regeneracyjnych poszczególnych kół funkcjonują- 
cych w zestawach kołowych. 


Pracę wykonano w ramach Badań Statutowych prowadzonych w Akademii 


Górniczo-Hutniczej w Krakowie na Wydziale Inżynierii Mechanicznej i Robotyki, 
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